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We conduct a stratigraphic paleoecological investigation at a Hudson River National Estuarine Research 
Reserve (HRNERR) site, Tivoli Bays, spanning the past 1100 years. Marsh sediment cores were analyzed 
for ecosystem changes using multiple proxies, including pollen, spores, macrofossils, charcoal, sediment 
bulk chemistry, and stable carbon and nitrogen isotopes. The results reveal climatic shifts such as the 
warm and dry Medieval Warm Period (MWP) followed by the cooler Little Ice Age (LIA), along with 
significant anthropogenic influence on the watershed ecosystem. A five-fold expansion of invasive 
species, including Typha angustifolia and Phragmites australis, is documented along with marked changes 
in sediment composition and nutrient input. During the last century, a ten-fold sedimentation rate 
increase due to land-use changes is observed. The large magnitude of shifts in vegetation, sedimentation, 
and nutrients during the last few centuries suggest that human activities have made the greatest impact 
to the marshes of the Hudson Estuary during the last millennium. Climate variability and ecosystem 
changes similar to those observed at other marshes in northeastern and mid-Atlantic estuaries, attest to 
the widespread regional signature recorded at Tivoli Bays. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Projected warming, drought, sea level rise, and salt intrusion 
(Bindoff et al., 2007; Christensen et al., 2007) make Atlantic estu- 
arine wetlands particularly vulnerable to future climate change. 
Alteration of wetlands within the estuary is especially of concern as 
wetlands serve many important ecosystem functions including 
hurricane protection, fish nursery habitat, and carbon storage. 
Understanding their resilience to past climatic shifts is paramount 
in providing a framework for their preservation and optimal 
function. Documenting ecosystem history provides multiple case 
studies of past ecosystem dynamics under quite varied environ- 
mental conditions (Jackson and Hobbs, 2009). However, such 
paleoecological data sets are rare, especially for the Hudson River 
Valley, New York, an extensive and intensively modified estuary. 
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The Hudson River watershed is an important resource (Findlay 
et al., 2002b) for over 10 million people in large population 
centers such as the metropolitan area of New York City as well as in 
small communities along the river (Burns et al„ 2007). Its rich 
natural ecosystem provides transportation and recreational func- 
tion as well as a vital water supply. In the past few centuries, the 
Hudson watershed has been significantly altered. Anthropogenic 
impacts, such as wetland destruction, invasive species expansion, 
pollution, and erosion are all of concern (Mills et al., 1996; Templer 
et al., 1998; Caraco et al„ 2000; Findlay et al., 2002a). 

While lower Hudson wetlands have been the target of recent 
paleo-investigation (Pederson et al., 2005; Peteet et al., 2006), the 
mid-Hudson region contains freshwater peatland archives that 
have not been investigated. Hence, there is a need to identify this 
base-line information to assess past anthropogenic activities and 
climatic patterns in relation to projected shifts in climate and 
vegetation in the Mid-Hudson Valley region. Important questions 
arise -1 ) How much has climate change and human impact altered 
these important wetland riparian ecosystems? 2) What shifts in 
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vegetation composition, sedimentation, and chemistry are 
apparent? and 3) Are there any detectable ecosystem and nutrient 
dynamic shifts after European settlement? 

In this study, we examine how climate and human impacts have 
influenced plant ecology, invasive species expansion, habitat loss, 
carbon storage, and nutrient dynamics over the past millennium 
based on the multiproxy analysis of sediment cores using paly- 
nology, macrofossil, sedimentological, and geochemical analyses at 
the National Estuarine Research Reserve (NERR) Tivoli Bays, Hud- 
son Estuary, NY. Comparison of these parameters from the marsh 
before and after European settlement allows identification of the 
recent human impacts on this estuarine ecosystem and provides 
insight regarding regional environmental changes. 

2. Study location 

The freshwater Tivoli Bays wetlands, 11 km north of the town of 
Red Hook, NY are the largest (6.88 km 2 ) Hudson River National 
Estuarine Research Reserve (HRNERR) site on the Hudson (Fig. 1). 
The site consists of a mixture of tidal marsh, tidal swamps, 
mudflats, and subtidal shallow habitats, hosting narrowleaf cattail 
( Typhci angustifolia), arrowhead ( Sagittaria spp.), jewelweed 
( Impatiens capensis), spatterdock ( Nuphar lutea ), water celery 
( Vallisneria americana), non-native purple loosestrife ( Lythrum 
salicaria), and non-native water chestnut ( Trapa natans ) along with 
a wide range of important fauna (Kiviat, 1983; Kiviat and Beecher, 
1991). Adjacent rocky hillside forests are composed of oak (Quer- 
cus rubrum, Quercus prinus), maple (Acer rubrum, Acer saccharum ), 
birch (Betula spp.) and ash (Fraxinus americana) with understory 
shrubs such as dogwood ( Cornus florida) and river margins of alder 


(Alnus sp.), willow ( Salix sp.). The upland bedrock is composed of 
marine sedimentary rocks (shale and graywacke) and quartz 
(Isachsen et al., 2000). A continental climate affects the Hudson 
Valley region with a small river maritime influence. The area 
receives average annual precipitation of 1280 mm, with average 
January and July temperature of -4.2 °C and 22.2 °C, respectively 
(NCDC, 2011) at nearby Poughkeepsie. 

3. Methods 

3.3. Sampling 

We analyzed four sediment cores from three locations (Fig. 1). 
The area north of Cruger Island Road (Fig. 1) is called Tivoli North 
Bay, while the area south of the road is called Tivoli South Bay. Site 
Tivoli A (Fig. 1) is 15 m north of Cruger Island Road (42° 1.85'N, 73° 
55.30'W). A woody tidal swamp today, it retains a mixture of 
arrowhead ( Sagittaria sp.), buttonbush ( Cephalanthus occidentalis), 
alder ( Alnus sp.), ash ( Fraxinus sp.), and maple (Acer sp.). Core 
05TivoliNBl (Fig. 1) was taken in August 2005 using a Hiller peat 
corer and was used for the top 50 cm retrieval to prevent peat 
compaction. 05TivoliNB2, a 1-m core, was taken in November 2005 
within 1 m of the first sampling site and retrieved using the 
modified Livingstone corer (Wright et al., 1984). These two cores 
are used for all analyses. 07Tivoli3, a core from Tivoli B, was taken 
with a Dashnowski corer from the second sampling site at the edge 
of the marsh in Tivoli North Bay (42° 1.98'N, 73°55.5'W), 12 m north 
of Cruger Island Road and 3 m from the railroad track (Fig. 1 ). This 
site is dominated by Typha angustifolia. The fourth core (Fig. 1), 
07Tivoli4 (42° 1.81'N, 73°55.59'W), was taken from Tivoli C, which 
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Fig. 1 . a) Hudson estuary watershed map, showing Tivoli Bays, mid-Hudson River, New York State with three other Hudson River National Estuarine Research Reserves sites, and 
b) Three Tivoli Bays sampling sites. The sites are located 11 miles north of town of Red Hook north and south of Cruger Island Road. 
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is 5 m west of the railroad in the South Bay and 20 m south of 
Cruger Island Road in T. angustifolia and Sagittaria spp.-dominated 
marsh. The cores were sub-sampled contiguously at 2 cm intervals, 
with an archive quarter stored at the LDEO repository. Approxi- 
mately 10 cc of homogenized samples from 05TivoliNBl and 
05TivoliNB2 were preserved wet for pollen and macrofossil anal- 
yses. The remainder of each sample was prepared for further 
analysis by drying at 105 °C and pulverizing with a mortar and 
pestle. 

3.2. Lithology 

Sediment lithology was analyzed based on visual inspection of 
color, texture, grain size, grain type, and fiber content using the 
USDA classification system. Sediment was delineated into different 
lithology types where sediment type abruptly changes based on 
these parameters. 

3.3. Loss-on-ignition (LOI) 

The dried samples were weighed and placed in a muffle furnace 
at 375 °C for 8 h (modified from Dean (1974).) The LOI is equal to 
sample weight lost as a result of combustion relative to the initial 
dry sample weight and representative of the amount of organic 
matter. 

3.4. X-ray Fluorescence spectroscopy 

Approximately 3 g of dried sediment was analyzed for chemical 
composition with X-ray Fluorescence spectroscopy (XRF) using 
Innov-X Alpha series 4000 XRF (Innov-X Systems, Woburn, MA) 
following protocols detailed by Kenna et al. (2011). Each analysis 
included two 120 s measurements using the soil protocol. Each 
sample was first analyzed in standard mode (Pb and Zn in this 
study) and then in Light Element analysis mode (K for this study). 

3.5. Gamma ray spectroscopy 

After being sealed for ~3 weeks to allow for short-lived 238 U- 
series radionuclides to re-establish equilibrium, approximately 
3—6 g of dried sediment from the top 50 cm of core 05TivoliNBl 
were analyzed at 2-cm intervals for gamma emitting radionuclides. 
The gamma counting system was comprised of a Canberra HpGe 
well detector model GCW2022 connected to a multi-channel buffer 
(Ortec model 919 Spectrum Master; 16,384 channels) and 
a personal computer running the MAESTRO® multi-channel 
analyzer emulation software package (Ortec). Detector efficiencies 
were determined by counting an International Atomic Energy 
Agency standard (IAEA-375) and U.S. Environmental Protection 
Agency standard pitchblende ore (SP-1). Corrections to the 
measured radionuclide counts included background, counting 
efficiency, geometry, branching ratio, and dry sample mass to 
obtain activity per dry weight gram. 210 Pb activities were corrected 
for self-attenuation (Appleby et al., 1992). Excess 210 Pb was calcu- 
lated by subtracting the 226 Ra from the total 210 Pb. Excess 210 Pb and 
137 Cs activities were then decay corrected to core collection date. 
The precision of the measurement was ±2—10%, and the back- 
ground levels were all <0.05 cpm. 

3.6. Radiocarbon dates 

Identified macrofossils were selected and radiocarbon dated 
using an Accelerator Mass Spectrometer (AMS) at Lawrence Liver- 
more National Laboratory. Radiocarbon dates were calibrated to 


calendar dates using the CALIB program, version 6.0 of Stuiver and 
Reimer (1993). 

3.7. Pollen, spore, and charcoal analysis 

Wet sediments were processed for pollen analysis with the 
standard HC1, HF, and acetolysis treatment, along with 7 and 150 pm 
filtration and centrifugation to extract the pollen and spore fraction 
(Faegri and Iversen, 1975). Exotic Lycopodium spores were added to 
calculate pollen concentration. Identification was based on pollen 
identification keys and LDEO reference pollen (Faegri and Iversen, 
1975; Kapp, 2000; Finkelstein, 2003). Typha angustifolia was care- 
fully separated from other Typha species and from Sparganium based 
on its shape and structure (Finkelstein, 2003). Phragmites australis 
pollen was separated from other Poaceae based on its smaller size 
than other Poaceae spp. in tidal marshes (Clark and Patterson, 1985). 
A minimum of 300 identifiable pollen grains for each sample were 
counted at 400x magnification. Pollen percentage was calculated 
based on total pollen sum, excluding spores. Charcoal counts 
included only charcoal pieces >250 pm 2 . Charcoal is presented as 
charcoal to total pollen (charcoakpollen) ratio. 

3.8. Macrofossil analysis 

Five cm 3 aliquots of wet samples were passed through a 500 pm 
sieve and rinsed with water. Material retained on the sieve was 
suspended in water and examined under a microscope at up to 80 x 
magnification. Macrofossils were identified using an extensive 
reference seed collection at LDEO and plant identification keys 
(Fernald, 1970; Martin and Barkley, 1973; Montgomery, 1977). 

3.9. Total C and N and stable C and N isotope analysis 

Dried and ground samples were weighed into tin cups (1—2 mg 
per sample). Three replicates of each sample were analyzed on 
a CHNS/O analyzer (2400 Series II, Perkin— Elmer, Boston, MA, USA) 
with combustion and reduction temperature of 925 °C and 640 °C 
for C and N content respectively. Acetanilide (1—2 mg) was used as 
standard and was run every 10 samples to ensure the machine’s 
stability. Carbon and nitrogen fluxes are calculated as %C (or 
N) x bulk density x sedimentation rate. 

Three other replicates of each sample were analyzed using a PDZ 
Europa 20-20 mass spectrometer with an ANCA-SL combustion 
system (Crewe, UK). The system measures total C, total N, 5 13 C and 
5 1S N. The delta (5) notation refers to the following ratio (R) of rare 
versus abundant isotope of either C or N with respective standard 
(PD-Belemnite and air respectively) according to Craig (1957), 
where 

- 1^1000 

A student t-test (two-tail) was performed to determine if two 
groups of dataset are significantly different. 

4. Results 

4.1. Lithology, LOI and sediment chemistry 

Four sediment types comprise the upper 100 cm of the cores 
(Fig. 2). From 100 to 66 cm, the sediment is composed of fine black 
peat muck with fibrous material. From 66 to 52 cm, the dark 
sediment is very fine. The lithogenic fraction of the top two zones is 
composed of well -sorted gray silty clay, though the lower zone is 
darker. A significant color change from dark gray to a lighter gray 
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Fig. 2. The lithology and chemistry of Tivoli Bays sediment cores, a) organic (black) and inorganic (gray) content in samples from Tivoli A, b) total K and c) Pb from XRF analysis of 
Tivoli A (circle), Tivoli B( triangle), and Tivoli C (diamond) follow inorganic matter trend with an increase after inorganic content increases (red line). (For interpretation of the 
references to colour in this figure legend, the reader is referred to the web version of this article). 


color occurs at 32 cm. Organic matter fluctuates between 0.8 and 
1.1 g cm 3 . Inorganic matter content increases abruptly at 52 cm 
(Fig. 2), and remains elevated to the surface. This corresponds to a K 
concentration increase above 50—60 cm at all three locations 
(Fig. 2), paralleling a major sedimentation change with higher 
inorganic matter input. Other lithogenic indicators, such as Ti (not 
shown), follow a similar trend as K. The concentration of trace 
metals, including total Pb (Fig. 2), of all sites also increases after 
50—60 cm depth, indicating that the lithogenic and chemistry 
change above these depths happened in the early 20th century 
(Nitsche et al., 2010). 

4.2. Chronology 

In order to estimate the sediment accumulation rates in the 
upper portion of the core, we used the constant flux: constant 
sedimentation (CFCS) excess 210 Pb model (Krishnaswami et al., 
1971; Robbins et al., 1978), fitting an exponential curve to the 
excess 210 Pb activity versus depth (Fig. 3a). Excess 210 Pb ( 210 Pb xs ) 
activities ranged between below detection and 3000 pCi/kg 210 Pb xs 
is complicated in this core. Evidence of disturbance in the profile 
can be observed at 20 cm. Values fall below detection at 32 cm, 
indicating perhaps a second discontinuity. Sedimentation rates 
estimated for intervals 0—20 and 20—32 independently are similar 
(0.8 and 0.74 cm/yr, respectively, Fig. 3a). 

The excess 210 Pb age model yields deposition ages between 1961 
and 1972 for intervals between 24 and 30 cm (Fig. 3b), which 
contain the 137 Cs maximum. The timing is consistent with global 
fallout deposition related to atomic weapons tests (Hardy, 1977; 
Monetti, 1996). We also observe above background levels of Pb and 
Zn (Figs. 2and 3), metals that are associated with anthropogenic 
industrial activities above 50 cm at all sites. The agreement 


between independent chronostratigraphic indicators provides 
a measure of confidence in the chronology above 32 cm and also 
suggests that the disturbance observed at 20 cm is minor. 

Deposition ages between 32 and 52 cm are less certain based on 
210 Pb. However, the continued presence of total Pb and total Zn 
(Fig. 3) allows us to constrain the timing to the early 20th century 
(Nitsche et al., 2010). Extrapolation of the excess 210 Pb derived 
sedimentation rate of 0.74 cm/yr yields a deposition age of ~ 1930 
for the discontinuity at 52 cm. 

The depth to age relationship below 52 cm is based on two 14 C 
dates for macrofossils collected at 85 and 99 cm (Table 1) which 
yields a sedimentation rate of ~0.05 cm/yr (Fig. 4). If one applies 
the 14 C based sedimentation rate, the deposition age of the base of 
the discontinuity at 52 cm is ~ 1740 AD (Fig. 4), which is consistent 
with the timing of the European Settlement of the early 18th 
century. Although there is evidence of some disturbance/disconti- 
nuity in the sedimentary record, information obtained from the 
application of multiple independent chronostratigraphic indicators 
is consistent, suggesting that sediments obtained from Tivoli Bays 
can provide a robust comparison of two distinct time frames - 
before and after the European settlement. The consistent deposi- 
tional patterns at three sites suggest that the lithogenic changes are 
widespread in the marsh. 

4.3. Pollen and spore assemblages and charcoal 

Based on visual inspection in Fig.5, the pollen stratigraphy from 
100 to 0 cm can be divided into four pollen assemblage zones. The 
pollen concentration profile yields a similar pattern and thus the 
concentration profile is omitted. Fig. 6 illustrates a simplified 
diagram of non-invasive, invasive, and woody wetland species to 
illustrate ecological compositional shifts. 
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Tivoli A Excess Pb-210 Profile Cs-137 Activity (Bq/Kg) 



Fig. 3. a) Excess 210 Pb and b) 137 Cs activity (blue circle), total lead (Pb, triangle) and zinc (Zn, cross). 210 Pb activity is plotted on a log scale and the profile follows a linear trend with 
discontinuity at 20 and 32 cm depths. The count abruptly drops below detection (B.D.) at 34 cm 137 Cs rises above background level at 40 cm and reaches a maximum value at 28 cm 
depth. Total Pb and Zn (b) based on X-ray Fluorescence analysis rise from the background above 60 cm and reach their maximum at 28 cm. (For interpretation of the references to 
colour in this figure legend, the reader is referred to the web version of this article.) 


4.3. 2. Zone 1: Quercus-Pinus-Betula-Carya-Osmunda (WO— 78 cm) 

Quercus (22—30%), Pinus (13—20%), Betula (7—10%), Carya 
(3—5%), and Osmunda (18—25%) reach maximum percentage in this 
zone. Poaceae (10—17.5%), Cyperaceae (12—22%) and Polypodiaceae 
(10—22%) percentages are also relatively high. Picea (2.5% average) 
is consistently visible in this zone. Charcoal increases abruptly in 
the upper half of the zone from 92 to 78 cm and is close to the 
maximum for the entire core. 

4.3.2. Zone 2: Quercus-Pinus-Tsuga-Cyperaceae-Poaceae- 
Polypodiaceae-Osmunda (78—52 cm) 

Quercus (18— 22%) and Pinus (10— 15%) remain the two dominant 
arboreal species. Populus reaches its maximum (8%) at 70 cm and 
Tsuga reaches maximum values. Pollen of herbaceous taxa account 
for 35—45%, with Cyperaceae (25—35%) increasing abruptly and 
Poaceae declining toward the top of the zone. Polypodiaceae 
reaches maximum values, and Osmunda remains high but declines 
toward the upper portion of the zone. Nymphaea was present 
throughout the zone. Charcoahpollen reaches the maximum value 
in the lower part of the zone at 74—72 cm. 

4.3.3. Zone 3: Quercus-Alnus-Poaceae-Ambrosia-Typha 
(52-36 cm) 

The arboreal fraction is comprised mostly of Quercus (10—15%) 
and Alnus (20—25%). Percentages of Tsuga and Pinus decline to 
minimum values of 5% in this zone. Betula declines at the 


Table 1 

AMS radiocarbon dates of Tivoli bays core. 


Lab 

Number 

Depth 

Materials 

dated 

Uncorrected 
14 C Age 3 

95.4% 

(2a) cal 
age ranges 
(area) 

Selected 
calendar 
age (AD) 

N74940 

84-86 

cm 

2 Rosaceae 
seeds and 1 
Carex seed 

965 ± 35 

1015-1160 

(85%) 

1087.5 

N85635 

98-100 

1 Carex seed 

1250 ± 90 

645-977 

(95%) 

811 


a With two standard deviation age range as calibrated by CALIB 6.0. 


beginning of the zone, but starts to increase at 40 cm. Ambrosia, 
Chenopodiaceae, Poaceae and Impatiens markedly increase, 
concurrent with a rise in Salix and Cephalanthus, while 
Polypodioceae and Osmunda sharply drop. 

4.3.4. Zone 4: Quercus-Pinus-Betula-Typha-Poaceae (36—0 cm ) 
Quercus and Pinus recover, with a rising trend reaching the 
maximum of 25 and 15% respectively. Betula (8%) recovers its high 
percentage as Ambrosia declines. Typha angustifolia (30—60%) 
becomes dominant. Salix, Fraxinus, Cephalanthus, Vitis, and Viola- 
ceae continue to increase. Lythraceae and Myriophyllum can 
consistently be observed throughout the zone. Other aquatic taxa, 
including Nymphaea and Potamogeton become noticeable. The 
charcoahpollen ratio declines to less than 30, but then slightly 
increases mid-zone. 


Calendar Year (AD) 

600 800 1000 1200 1400 1600 1800 2000 2200 



Fig. 4. Age-depth relationship of Tivoli A. An extrapolation based on AMS 14 C age 
below 52 cm of the core at 0.05 cm/yr suggests that the 52 cm boundaiy is dated back 
to 1740 AD. An extrapolation using 0.74 m/yr sedimentation rate in the upper 52 cm 
suggests that the 52 cm depth is dated to 1930 AD. 


S. Sritrairat et al. / Estuarine, Coastal and Shelf Science 102-103 (2012) 24—35 


29 



Zone 


4 


3 


2 


1 


Light gray silty clay 


Dark gray silty clay 



Black fine muck 



Black peat muck 


Fig. 5. Tivoli Bays pollen percentage profile. Percentage is calculated in relation to total pollen sum and is categorized by vegetation type. 


4.4. Macrofossil assemblages 

The macrofossil diagram is presented in Fig. 7 utilizing the same 
zones as in Fig. 5. Very few macrofossils are found in the lowermost 
zone (Fig. 7). Carex cf. aquatilis seeds throughout zone 2 are 
concurrent with Cyperaceae pollen increases in zone 3 (Fig. 5). 
Seeds of weedy species, including Asteraceae, Poaceae, Polygonum, 
Impatiens capensis and Chenopodiaceae are found primarily in the 
top zones (zone 3 and 4). Flowever, both pollen and seeds of I. 
capensis were observed in earlier zones as well. 

The macrofossils provide the advantage of confirming the 
presence of some plants as well as aiding identification at the 
species level. For example, Cornus obliqua seeds (zone 3) were 
found, though their pollen were not recognized in the pollen 
profile. Alnus bracts were found in zone 1, 2 and 4, even though 
Alnus pollen percentage is the highest in zone 4 only. Violet seeds 
are also found at several depths throughout, even though Violaceae 
pollen (Fig. 5) is the highest in zone 4. 


4.5. XC, XN, and C/N ratio 

Hydrogen to carbon (H/C) atomic ratio ranges from 0.98 to 1.45, 
indicating a low inorganic carbon fraction. Percent carbon and C/N 
ratio both show overall declines upward (Fig. 8). In zone 1, percent 
carbon decreases from 39% to 28% in this lowermost part of the core 
(Fig. 8). The C/N ratio in this zone from 98 to 80 cm is significantly 
higher than the rest of the core (P < 0.001 ), with a mean value of 
20.6 ± 1.5. Percent carbon declines slightly in zone 2 (26.6 ± 1.1%), 
while percent nitrogen (1.66 ± 0.13%) and the C/N ratio 
(16.0 ± 0.7%) are very stable during this interval. In zone 3, percent 
carbon gradually drops by half, from 25% to 12% over this interval. 
Percent nitrogen follows a similar decreasing trend, changing from 
1.6% to 0.7%. There is a slight decline of the C/N ratio where charcoal 
also drops, but both increase at the close of the interval. In zone 4, 
percent carbon content fluctuates around 11.0 ± 1.9%. Percent 
nitrogen is highly variable (0.668 ± 0.107%), while the C/N ratio is 
stable (16.5 ± 0.7%). 


Calendar Year (AD) 


% Pollen and Spores 


2005 


1966 ± 5 


1087 ± 72 
811 ± 166 


0 20 40 60 80 100 



Fig. 6. A comparison of the sum of percent pollen and spores of three plant groups. 1) non-invasive wetlands species (open circle, e.g. Cyperaceae spp., Osmunda sp, and Poly- 
podiaceae); 2) invasive (black circle, e.g. Phragmites australis, Lythraceae, and Typha angustifolia ); 3) woody wetland species (triangle, e.g. Salix spp. and Vitis spp.) at Tivoli Bays. 
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Light gray silty clay Dark gray silty clay Black fine muck Black peat muck 


Fig. 7. Macrofossils identified in Tivoli North Bay core. The fossils are larger than 150 pm and are seeds unless specified otherwise, in the unit of count per 5 cm 3 . 


4.6. Stable C and N Isotopic analysis 

Stable isotope C and N show significant fluctuations throughout 
the core (Fig. 8). The 8 13 C is stable at lower depths from 93 to 57 cm, 
with the mean value of -26.52 ± 0.37%,, Above 57 cm depth, the 
5 13 C increases, until reaching a maximum value of § 13 C 
(-22.82 ± 0.46% o ) at 47 cm. From the maximum, 5 13 C declines 
toward the top of the core to the minimum of -29.02 ± 0.11% o for 
all depths above 93 cm. 

From 100 to 78 cm, 5 1S N (Fig. 8) fluctuates around a low mean of 
1.47 ± 0.67%„ then drops to 0.40 ± 0.30% o at 78 cm. From 78 cm up 
to 38 cm, 5 15 N increases. From 38 cm upward to 12 cm, 5 15 N is 


stable and relatively enriched (4.06 ± 0.38% o ). However, it drops to 
3.5%„ around 9—5 cm before returning to the most enriched values 
at the surface (4.7% c ). 

5. Discussion 

5.1. Multiproxy interpretation of Tivoli Bays environmental change 

The mid-Hudson Valley marshes, as evidenced by this study, 
reflect enormous impact from humans. Vegetation composition, 
invasive species expansion, nutrient dynamics, and sedimentation 
patterns significantly shift at the boundary of the European 



Light gray silty clay Dark gray silty clay Black fine muck Black peat muck 


Fig. 8. Carbon and nitrogen parameters at Tivoli A. The graph shows %C, %N, %H from a CHNS/O analyzer, calculated C and N fluxes, and 5 13 C and 8 15 N from a spectrometer of bulk 
sediment core. 
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settlement at 52 cm (~1700 AD) (Fig. 7). Pollen percentage of 
native wetland species, including ferns and sedges, decreases by up 
to 75% from the pre-European settlement period toward the 
present. Total percent of invasive species, such as Typha angustifolia, 
Phragmites australis, and Lythrum salicaria (Fig. 6) increases five- 
fold. Invasive species expansion is clearly anthropogenic in most 
cases, as many species were not present or were present at low 
percentage prior to European settlement. Weedy and woody 
species also increase by three times their pre-settlement percent- 
ages, as supported by disturbance and drier conditions from a 10- 
fold (Fig. 4) higher sedimentation rate and higher inorganic matter 
input (Fig. 2) during the last century. Based on sediment chemistry, 
this inorganic matter increase is uniform throughout the marsh, as 
Ti and I< increase at all sites when total Pb and Znare above back- 
ground level (Figs. 2and 3), indicating influences from humans. 
Nutrient and carbon concentration, fluxes, and isotopes signifi- 
cantly shift (Fig. 8). Using the vegetation zones previously estab- 
lished based on pollen assemblages, detailed discussion about each 
time zone follows. 

5.1.1. Pre-European settlement with fires (A.D. ~826—1310) 

Forest composition at about 1000 AD is dominated by Quercus 
and Pinus (Fig. 5). Species that have an advantage in disturbed areas 
such as Betula and Poaceae are likely indicative of land disturbance 
from fires, as higher charcoal is observed. High percentage of Carya, 
a warmth-loving species (Fowells, 1965), supports an increase in 
temperature. The depth of 84 cm that is dated to 1087 =t 72 AD 
corresponds to the lower charcoal maximum, a feature that is also 
found in other Hudson River marsh cores at Piermont (Pederson 
et al., 2005) and Iona (Peteet et al., 2006) and represent the 
warm, dry Medieval Warm Period (MWP.) This high charcoal peak 
is likely a result of a regional Hudson Valley MWP recorded on 
a larger spatial scale in other parts of North America and the globe. 
For example, Chesapeake records suggest drier conditions prior to 
the Little Ice Age (LIA) (Brush, 1986, 2001 ; Cronin et al., 2000, 2010; 
Cronin and Vann, 2003) and the western US is known to have 
suffered severe droughts based on tree ring records (Cook et al., 
2004). While there is a debate if the MWP is a global phenom- 
enon as the warming is not synchronized at all sites around the 
globe, many paleoclimatic records suggest widespread climatic 
anomalies, such as parts of Europe (Mangini et al., 2005), Tasmania 
(Cook et al., 1991), Asia (Yang et al., 2002), and Africa (Alin and 
Cohen, 2003) during the same time period. 

During this time period prior to the European settlement, the 
sedimentation rate of 0.05—0.06 cm/yr is much lower than the 
sedimentation rate at the brackish downriver Piermont marsh 
(Fig. 1, 0.3 cm/yr (Pederson et al., 2005)) and is lower than most 
freshwater tidal marshes (Neubauer, 2008). The sediment is highly 
organic throughout, suggesting low contribution of any upland 
sediment source. 

If local drought and fire occurred during this period, enhanced 
erosion from the uplands would likely increase lithogenic input. 
This may explain the slight decline in percent carbon (Figs. 2 and 8), 
but the complexity of plant type in net carbon sequestration makes 
the connection obscure. An increase in the cold indicator, Picea 
pollen, overlapping the charcoal maximum (Fig. 5) may be a result 
of increased erosion from the sites further north. Despite the higher 
charcoal counts compared with the modern deforestation period, 
the forest and marsh composition of this pre-industrialization 
interval does not show as much disturbance in species composition 
as during the subsequent European settlement and industrializa- 
tion. The presence of aquatic species such as Nymphaea during this 
presumably warm and dry period suggests that the local plant 
community may not have been significantly affected by fires. 
However, the rise of ferns ( Osmunda and Polypodiaceae) and 


Cyperaceae may indicate a drop in water table at the site that allows 
the colonization of these genera near shore on exposed riverine 
mud. Fern expansion was observed at both Piermont and Tivoli 
supporting a regional hydrological shift. 

The 5 13 C signature in marsh sediments reflects the organic 
matter source, decomposition, and environmental factors (Chmura 
et al., 1987; Chmura and Aharon, 1995). Chmura et al. (1987) found 
that 8 13 C in sedimentary carbon in Louisiana marshes was mainly 
derived from the local dominant vegetation. Since some plant 
groups have a distinct 8 13 C range and because marsh vegetation has 
significantly changed at Tivoli based on pollen/macrofossil stra- 
tigraphy, vegetational composition is likely the most significant 
cause of 8 13 C fluctuations at this study site. There exist both C3 and 
C4 Cyperaceae which are the dominant pre-European marsh 
vegetation in this zone. Elevated 8 13 C signature observed in zone 1 
(Fig. 8) suggests possible expansion of C4 Cyperaceae as C4 
Cyperaceae are reported to be veiy enriched in 8 13 C ( — 12.1 to 
-15.9% 0 ) (Chmura et al., 1987; Chmura and Aharon, 1995). 

The drier condition of the wetland can result in reduced frac- 
tionation and decreases in S 13 C from vegetation change, while the 
decomposition may increase the 13 C content (Chmura et al., 1987; 
Malamud-Roam and Ingram, 2004). While the S 13 C is more 
enriched than in other zones, supporting high decomposition in 
a dry condition, the inferred lower water level during the warm and 
dry MWP at Tivoli does not seem to be evident based on S 13 C as the 
value is very stable and does not track the charcoal count. Tivoli 
remained fresh, and the vegetation type in the marsh may not have 
been significantly altered in comparison to subsequent vegetational 
change. 

There is a higher fluctuation of 8 15 N around lower mean value of 
1.47 =t 0.67% o below 80 cm (~1000 AD) (Fig. 8) when charcoal is 
high (Fig. 5). Burning might have caused the N fractionation and 
depleted 15 N (Ansley et al., 2006) during charcoal maxima. Drier 
conditions would also allow an increase in erosional terrestrial 
nitrogen input and higher nitrification rate. Both processes can 
contribute to a heavier 15 N signature than the products of denitri- 
fication which occur in conditions of greater anoxia in the sediment 
(Brandes and Devol, 2002). 

5.3.2. Pre-European forest and climate transition (A.D. 

~ 1310 1700) 

Quercus and Pinus continue to dominate the forest in this region 
before European settlement (Fig. 5). Quercus bicolor which is found 
in the woody Tivoli swamp today might have contributed to total 
Quercus pollen in addition to an upland signal (Kiviat, 1983). Pinus 
strobus is likely the dominant Pinus species, based on modern 
assemblages. 

Kiviat and Beecher (1991) estimate from an old 1850 map that 
North Bay was probably dominated by subtidal communities and 
deep open water. This is supported by our pollen stratigraphy that 
there were high amounts of aquatic pollen and spores, including 
Nymphaea and Potamogeton. Cyperaceae, Osmunda, and Poly- 
podiaceae could have resided around the marsh perimeter and on 
the deltas off several creeks that are present. 

Carex cf. aquatilis seeds (Fig. 7) confirm that Cyperaceae was an 
important marsh component during the pre-European settlement 
time. This sedge grows well in shallow pools, pond, and near stream 
margins (Fernald, 1970) and its presence suggests that Tivoli Bays 
may have been wetter and more open to river flow. Royal fern 
(Osmunda regalis ) is the single Osmunda species in Tivoli Bays today. 
However, Osmunda spores observed in this zone are comprised of 
0. regalis and at least one other species based on spore texture and 
size. Thus, higher Osmunda diversity before European settlement is 
suggested, including possibly O. cinnamomea and/or 0. claytoniana. 
In contrast to the locally moist marsh condition, a relatively high 
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percentage of Carya, Liquidambar, Populus , and charcoal count 
indicate that regional climate was probably drier as well as warmer 
at the beginning of this zone during the MWP transition. As char- 
coal dropped off, mesic species such as Tsuga began to rise, 
consistent with a transition from drier to wetter conditions. This 
wetter and cooler period during the L1A is also observed in Pier- 
mont (Pederson et al., 2005) and Chesapeake Bay (Brush, 2001; 
Cronin et al., 2010). 

Slightly lower C/N ratio early in this zone may be a result of 
preferential carbon loss from fire. 8 13 C continues to be stable and 
similar to the previous zone, with values higher than in other zones. 
Similar to zone 1, this elevated § 13 C supports the continuation of 
a high population of Cyperaceae with enriched 13 C. S 15 N rebounds 
and continue to increase as charcoal and fire occurrence decline. 

5.2.3. European settlement, industrialization, and landscape 
changes (A.D. ~1700—<1960s) 

Historic records indicate that the Europeans bought Cruger 
Island and Tivoli Bays from Native Americans in A.D.1681, but 
significant logging might not have started until 1688 (Klein, 2001; 
Philip, 2008). The increase in Ambrosia has been used to success- 
fully indicate the timing of European settlement in the Hudson 
Highlands and the Hudson Valley (Maenza-Gmelch, 1997 ; Pederson 
et al„ 2005). Thus, the Ambrosia rise at 54 cm (Fig. 5), if reflecting 
regional settlement, suggests the timing of the permanent Dutch 
settlement starting in the late 17th century (Philip, 2008) which 
agrees with the age estimation from radiocarbon dates. The major 
characteristic of this interval is a sharp decline of forest compo- 
nents Quercus, Pinus, Carya, Tsuga, and Liquidambar due to logging 
with settlement. Charcoal increases probably indicate human- 
induced fires. The expansion of Betula supports a disturbance 
scenario. Early-colonizing herbaceous species, including Poaceae, 
Chenopodiaceae, Impatiens, and Violaceae became much more 
abundant (Fig. 5) in response to land clearing. Higher percentage 
mesic and northern species such as Tsuga prior to this zone and 
Fagus in the upper part of the zone suggests a moist and cooler 
climate. Historic records indicate cold winters with high precipi- 
tation in the 1800s (Cook and Jacoby, 1977; Klein, 2001). 

Lithological change and the reduction of percent carbon are 
observed in this zone (Fig. 5). At all three Tivoli locations, higher 
inorganic matter input and higher K, which are indicative of 
increase of inorganic matter, support the shift in regional sedi- 
mentation pattern (Fig. 3). This sediment composition shift is likely 
due to higher erosion from deforestation during the European 
settlement and following industrialization. A similar trend is 
observed during this time frame at Piermont Marsh (Pederson 
et al., 2005), suggesting regional land-use changes in the estuary. 
As inorganic deposition rate was higher during this period, it 
increased the sedimentation rate and raised the elevation of the 
marsh, making the soils drier with improved drainage. A reduction 
of aquatic and wetland plants (Nymphaea, Potamogeton, Cyper- 
aceae, Osmunda, and Polypodiaceae) and replacement with 
shrubby wetland plants ( Alnus , Cornus, Cephalanthus and Vitis) 
supports this inference. Typha angustifolia, along with other inva- 
sive species, began to colonize the marsh (Figs. 5and 6). 

C/N ratios during this settlement interval are lower than in other 
zones (p < 0.04). This reduction indicates higher contribution from 
low C/N sources. Low C/N sources include degraded “old soil,” 
aquatic plants, and lower aquatic autotrophs such as plankton and 
algae. The deforestation could have reduced fresh terrestrial plant 
input. Subsequent erosion can bring in more well-oxidized upland 
soil, which has a lower C/N content (higher nutrient) than new 
organic material (Fisher et al., 2003). This initiation of nutrient 
increase at Tivoli possibly contributed to an increase in invasive 
species in the wetland. 


There is a significant decrease in § 13 C (p < 0.001 ) and increase of 
5 1S N (p < 0.01 ) compared to the rest of the core just after the 
European settlement at 48 cm (Fig. 8). The changes occur at the 
same depth as the lithological change and a major shift in pollen 
signifying disturbance. This concurrence supports the hypothesis 
that 5 13 C and 5 1S N shifts are due to anthropogenic activities. The 
depletion of 5 13 C might be a result from a shift in other sediment 
sources, such as possible increased input from erosion after land 
clearing or the increase in algae (minimum algae § 13 C of -17.9% 0 , 
(Chmura et al., 1987)) growth due to nutrient fertilization from 
upland runoff, which is also supported by the C/N ratio. 

The change and increase in different nitrogen sources is likely 
the main factor that drives the increasing 5 15 N trend. Human and 
animal waste input (enriched in 15 N (Aravena et al., 1993)) to 
estuarine systems can result in higher 8 1S N in marsh vegetation 
(Cole et al., 2004). Increasing population and livestock in the 
watershed after the European settlement may have increased 
human and animal waste input which results in the observed 
higher 15 N. 

5.2.4. Reforestation and marsh restoration (<1963— present) 

Sedimentation rate in this zone at 0.74 cm/yr is more than 10 
times the sedimentation rate prior to the European settlement, 
probably as a result of land-use changes and estuarine disturbances. 
The rate is higher than that of current sea level rise (0.31 ± 0.07 cm/ 
yr in 1993—2003) and it falls within the range observed sedimen- 
tation rates in fresh tidal marshes from various locations (Neubauer, 
2008). Regrowth and reforestation is indicated by the resurgence of 
Quercus, Pinus, and Betula concurrent with the Ambrosia decline in 
the second half of the twentieth century. However, forest compo- 
sition is different from pre-settlement as Pinus is less abundant and 
Betula expands. Vacant and open agricultural land may also result in 
the continuous growth of weedy species. 

In contrast to the Pinus and Quercus reforestation, Tsuga sharply 
declined and remains low up to the present. In contrast to our 
pollen record suggesting lack of hemlock recovery, at some New 
England sites an increase in Tsuga canadensis has been observed as 
a result of replanting and a decline in logging during recent decades 
(Orwig and Foster, 1998; McLachlan et al., 2000). 

Vegetation of the marsh itself also differs from pre-European 
time (Figs. 5—7). Woody species, including Salix, Fraxinus, and 
Cephalanthus, continue to increase, while Cyperaceae, Osmunda, 
and Polypodiaceae never reach their pre-settlement percentage. 
Typha angustifolia becomes one of the dominant marsh plants, with 
more than five times higher pollen percentage than during the pre- 
settlement when its pollen was scarce. Tivoli plant surveys in the 
1970s and 1980s indicate that T. angustifolia has been a major 
wetland species for several decades, while woody plants were 
expanding toward the river (Westad and Kiviat, 1986). 
T. angustifolia and T. x glauca have become prominent species in 
North American wetlands and are often considered as weedy, 
unwanted or invasive with negative impacts on the ecosystem as 
they spread and develop monocultures, reducing plant diversity 
(Grace and Harrison, 1986; Houlahan and Findlay, 2004; Shih and 
Finkelstein, 2008; Vaccaro et al., 2009; Rippke et al., 2010). Based 
on herbarium and pollen records, Shih and Finkelstein (2008) 
concluded that Typha spp. has a high invasive tendency in North 
America. Thirty five out of 45 sites North America showed an 
increase in T. angustifolia ( Typha/Sparganium monads pollen type) 
in the past 1000 years, concurring with the increasing trend of this 
species from this study. Piermont Marsh paleoecology indicates 
a similar expansion of T. angustifolia (Pederson et al„ 2005) in 
recent decades. 

Vitis labrusca is very common on woody plants, so Vitis increase 
may be a result of the expansion of woody species. The presence of 
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Asteraceae, Poaceae, Polygonum, and Chenopodiaceae seeds indi- 
cate land disturbance and invasion of weedy species which is also 
observed in the pollen and spore profile. Higher sedimentation rate 
likely plays a role, as the emergent area of North Bay has increased 
from 1934 to 1995 near the northwest boundary and along Cruger 
Island Road based on historical maps. 

The most abundant herbs at Tivoli during the last three decades 
are Impatiens capensis ( Impatiens biflora), Lythrum salicaria, Polyg- 
onum spp., Osmunda regalis, Pilea pumila, and Viola cucullata 
(Westad and Kiviat, 1986). In our record, there may be up to a four 
times reduction in Osmunda spp. while the percentage of Vitis 
pollen is at least 3 times higher than pre-settlement time, though 
the seed concentration is stable throughout the core. Small size 
Poaceae pollen, likely representing Phragmites australis, shows 
a similar trend as it is invading Tivoli (Kiviat, 1983). Small size 
Poaceae pollen was present before the settlement, but a minor 
component. This is supported by other studies suggesting that 
P. australis was native to North America for at least 3000 years 
(Orson et al., 1987) based upon rhizome identification. Its aggres- 
sive invasion over the last 150 years in the Mid-Atlantic region 
(Galatowitsch et al., 1999; Rice et al., 2000; Silliman and Bertness, 
2004) is attributed to of the introduction of a more aggressive 
hybrid strain from Eurasia during the early 19th century 
(Saltonstali, 2002). By DNA sequencing of herbarium collections, 
the imprint of the introduced strain was shown to be present in the 
Mid-Atlantic region before 1910. By 1940, only the hybrid strain has 
been identified (Saltonstali, 2002, 2003). Land disturbances and 
nutrient enrichment have also been linked to the expansion of this 
species (Rice et al., 2000; Silliman and Bertness, 2004). 

The invasive Lythrum salicaria has been reported to cover large 
areas during the past three decades (Kiviat, 1983; Westad and 
Kiviat, 1986; Kiviat and Beecher, 1991). As seen from the pollen 
record, Lythrum was not significantly present until after European 
settlement. Myriophyllum pollen, consistently found at a low 
percentage in the upper zone, is also likely the introduced 
M. spicatum, the only Myriophyllum species in Tivoli (Westad and 
Kiviat, 1986). Impatiens pollen and Impatiens capensis seeds 
increase after the European settlement in similar nature to the 
increase of other invasive species. The expansion of this species in 
the fossil record appears unique to Tivoli. This species is native to 
North America and is considered to be non-invasive in the 
northeast United States, but it is invasive in Europe (Tabak and 
Wettberg, 2008). Hybridization of the species with other non- 
native Impatiens has not been documented in Northeastern 
United States, but Tobak and Wettberg suggest there might be 
a threat of such hybridization. The increase of the species might 
also be due to the change in environmental factors that favor its 
expansion and thus close attention to the spread of the species is 
warranted. 

Throughout the upper portion of the core, S 13 C linearly 
decreases to below the average pre-historic level. Compilation of 
the § 13 C of marsh vegetation reveals that modern invasive and 
weedy species (i.e. Typha x glauca (-27.7% 0 ), Phragmites australis 
(-24.6 to -29.4% 0 ), I. biflora (Impatiens capensis ) (-27.4% 0 ) are 
highly depleted in 13 C (Chmura and Aharon, 1995; Hornibrook et al„ 
2000) as are most upland trees. Because terrestrial upland and 
marsh species have more depleted 5 13 C than aquatic plants, 
decreasing 5 13 C is consistent with the observed woody species 
invasion into Tivoli Bays from the pollen and spore record (Fig.6). 
Degradation of organic matter can result in an enrichment of 8 13 C in 
the sediment, though selective decomposition might lower 5 13 C as 
more 13 C — depleted organic components such as lignin are better 
preserved (Chmura et al., 1987; Hornibrook et al., 2000). This can 
cause an additional effect of lowering the 8 13 C and the gradual 5 13 C 
decrease near the top of the core as labile-organic matter 


undergoes more rapid decomposition than old organic matter 
(Hornibrook et al., 2000). The 5 13 C trend somewhat resembles that 
of % organic matter profile (r 2 = 0.7, p < 0.01 ) in that 5 13 C is more 
negative with lower organic matter content as higher lithogenic 
input and sedimentation rate causes the invasion of woody species. 

5 1S N often parallels wastewater/run off input as well as eutro- 
phication increase in S 15 N in macrophytes and marsh sediment of 
the Hudson estuary and Chesapeake Bay (Cole et al., 2004, 2005; 
Abreu et al., 2006; Church et al., 2006). In this zone, S 15 N remains 
elevated in comparison to the pre-European settlement period, 
likely due to the continuation of increased wastewater input. In 
addition, nitrogen consumption in the United States has increased 
from near 0 tg/yr in the 1940’s (Vitousek et al., 1997), to 2.74 tg/yr in 
1960 which then exponentially increased to 13.2 tg/yr in 2007 
(Erseus, 2008). Significant increase of agricultural and sewage input 
to the Hudson Estuary during the last 50 years likely caused 
eutrophication in the estuary (Howarth et al., 1996). Nitrogen and 
carbon flux exhibits a significant increase in the last 50 years 
(Fig. 8), supporting this view. 

Anoxic conditions can increase denitrification, and thus enrich 
15 N in the sediment further. Evidence of anoxia is also seen as higher 
carbon and nitrogen fluxes (Fig. 8). Thus, the profile of higher 8 15 N 
after the European settlement at Tivoli supports increasing fertilizer 
usage, wastewater input, and eutrophication. Similar increase in the 
degree of eutrophication during the same time scale is observed in 
many other estuaries, including Chesapeake Bay (Cooper and Brush, 
1993), Narragansett Bay (Hubeny et al., 2009), and the Great Lakes 
(Rippke et al., 2010). The shift in 8 1S N begins at the same time that 
Typha angustifolia rises and may be linked. This correlation supports 
the hypothesis that T. angustifolia expansion is fostered by eutrophic 
and nutrient-rich conditions (Vaithiyanathan and Richardson, 1999; 
Tanaka et al., 2004). 

5.2. Tivoli Bays in a regional context 

To place the regional ecological changes of the Hudson marshes 
into a perspective based on historical records, Swaney et al. (2006) 
suggest that the Hudson watershed was probably heavily forested 
prior to 1609. Land along the river was colonized by 1700 AD 
because the land was easily accessible by boats (Swaney et al., 
2006). Swaney et al. also suggest that land clearing and logging of 
any region along the river took place early, which would explain 
pollen evidence for arboreal decline in the 17th— 18th century at 
Tivoli and other Hudson River sites (Pederson et al., 2005; Peteet 
et al., 2006). Palynological records from Hudson Highlands ponds 
(MaenzaGmelch, 1997) and nearby New Jersey lake (Toney et al., 
2003), suggest similar tree declines, and herbs and shrubs 
increase around 1700 AD. The tanning industry started in the 1700s 
and peaked in the mid 1800s, when the primary resource, Tsuga, 
declined. Both logging and tanning collapsed from exhaustion of 
resource and market by the mid 1850s. In 1880, 68% of the Hudson 
watershed was farmland, but it has been 70% reforested by the 
1990s (Swaney et al., 2006) and is evident as tree species recovery 
in pollen records at other sites (MaenzaGmelch, 1997; Pederson 
et al., 2005). This would explain the rise of tree species at all 
sites. As expected, pollen percentages of trees are lower than the 
pre-European values. 

Westad and Kiviat (1986) suggested that construction of the 
Cruger Island Bridge and railroad in 1835 and 1851 may have 
contributed to higher sedimentation rates in the bays. This devel- 
opment fragmented the bays and slowed down tidal flushing. From 
our record, a significantly higher inorganic content near the surface 
at all three Tivoli locations indicates a major landscape change. The 
railroad construction is the most likely cause of sediment trapping. 
Higher sedimentation rates may have caused the observed increase 
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in drier and woodier species. However, the higher percent of 
inorganic matter clearly increased when Ambrosia started to rise at 
Tivoli A, presumably in the eighteenth century. Therefore, it 
appears that higher inorganic input resulted from the increase in 
land erosion with settlement somewhat earlier, rather than solely 
from the subsequent local bridge, road, and railroad constructions. 

From this record, Tivoli marsh has stored a significant amount of 
organic matter, as the %C in the sediment ranges from 10 to 40%. 
However, the level drops off after the European settlement. In 
comparison, the downriver Piermont core has twice as much %C 
(20%) at the surface, and the salt marsh vegetational composition 
dominated by Spartina patens / Distichlis spicata there is quite 
different. The LOI at Piermont declined during settlement as the 
inorganic weight fraction first increased due to erosion but then 
recovered in the last century as reforestation led to erosional 
decline (Pederson et al., 2005). A consistently high inorganic frac- 
tion at Tivoli in recent centuries may be a result of continued local 
landscape change rather than a regional effect. 

Cook and Jacoby (1977) concluded from a tree ring study that 
drought during the 1960s was the most severe since 1730 in the 
Hudson Valley. However, charcoal concentration near the top of this 
core is five times smaller than during the 11th century maximum. 
Thus, in comparison, the MWP was probably much drier than recent 
centuries and additionally humans have suppressed fire. 

6. Conclusions 

In answer to our research questions, the Tivoli Bays paleorecord 
demonstrates an example of enormous anthropogenic impact of 
European settlement on an estuarine ecosystem at an NERR site. 
While upland forest composition recovers somewhat from the pre- 
European settlement period, wetland species composition, 
lithology and function do not recover their pre-European state. 
Sedimentation rate at Tivoli Bays during the last two centuries 
(0.74 cm/yr) is more than 13 times higher than the pre-settlement 
period (0.05 cm/yr), indicating significant impacts from human- 
induced land-use changes. A concurrent significant increase in 
woody species ( Salix spp., Vitis spp., and Alnus spp.) suggests a drier 
habitat from regional increase of inorganic sediment input. Local 
marsh vegetation at Tivoli Bays prior to European impact was 
dominated by sedge (Cyperaceae). Ferns such as Osmunda and 
Polypodiaceae were also well-represented. Species with high 
invasive tendency such as Phragmites australis, Typha angustifolia, 
Lythrum salicaria, and Impatiens capensis all increased up to 5 fold 
during the last 50—100 years and replaced pre-settlement sedges 
and ferns. The expansion of these species parallels pollen records in 
the Hudson Estuary and across North America, confirming the 
invasiveness of these species. Significant 5 13 C decline after the 
European settlement at Tivoli Bays confirms the vegetation 
compositional shifts. C and N fluxes, 8N, %N profiles support 
significant anthropogenic influences on nutrient dynamics of the 
wetland. A significant increase in wastewater input, fertilizer usage, 
and anoxia are likely linked to invasive species expansion. The 
finding confirms regional MWP with drought leading to fire and 
a cooler LIA climate. However, vegetation shifts as a result of these 
climatic events are smaller than those from human activities. 
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